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Abstract. We model the population of double white 
dwarfs in the Galaxy and find a better agreement with 
observations compared to earlier studies, due to two mod- 
ifications. The first is the treatment of the first phase of 
unstable mass transfer and the second the modelling of 
the cooling of the white dwarfs. 

A satisfactory agreement with observations of the local 
sample of white dwarfs is achieved if we assume that the 
initial binary fraction is ~50 % and that the lowest mass 
white dwarfs [M < 0.3Mq) cool faster than the most 
recently published cooling models predict. 

With this model we find a Galactic birth rate of close 
double white dwarfs of 0.05 yr -1 , a birth rate of AM 
CVn systems of 0.005 yr -1 , a merger rate of pairs with a 
combi ned mass exceeding the Chandrasckhar limit (which 



have experienced at least two phases of mass transfer and 
thus provide good tests for theories of binary evolution. 
Second it has been argued that type la supernovae arise 
from merging double CO white dwarfs ( Webbink 1984 ; 
[ben fc Tutukov 1984). Thirdly close double white dwarfs 



may be the most important contributors to the gravita- 
tional wave signal at low frequencies, probably even pro- 
ducing an unresolved noise burying many underlying sig- 
nals (Evans et al. 1987; Hils ct al. 1990). A fourth reason 
to study the population of double white dwarfs is that in 
combination with binary evolution theories, the recently 
developed detailed cooling models for (low-mass) white 
dwarfs can be tested. 

The formation of the population of double white 



J> . may b ; progenitors of SNe la) of 0.003 yr 1 and a for- 



mation 



dwarfs has been studied analytically by Iben & Tutukov 



rate of planetary nebulae of 1 yr . We estimate 



(1986a 
(1988|) 



1987) and numerically by Lipunov & Postnov 



Tutukov fc Yungelson (1993, 1994); Yungelson 



the tot al number of double white dwarfs in the Galaxy as 
2.5 x 10 8 . 

In an observable sample with a limiting magnitude 
Vn m = 15 we predict the presence of ~855 white dwarfs 
of which ~220 are close pairs. Of these 10 are double CO 
white dwarfs of which one has a combined mass exceeding 
the Chandrasekhar limit and will merge within a Hubble 
time. 



et al. (1994| ); [Han ct al. (1995| ); [Iben et al. (1997[ hereafter 
[TY97|), and [Han (1998|,Tereafter |HAN98[). Comparison 



Key v jords: stars: white dwarfs - stars: statistics - bina- 
ries: close - binaries: evolution 



1. Introduction 

Close double white dwarfsQ form an interesting popula- 
tion for a number of reasons. First they are binaries that 



between these sudies gives insight in the differences that 
exist between the assumptions made in different synthesis 
calculations. 

Following the discovery of the first close double white 
dwarf ( gaffer et al. 1988| ) , the observed sample of such sys- 
tems in which the mass of at least one component is mea- 
sured has increased to 14 ( JVIaxtcd fc Marsh 1999 ; Maxted 
ct al. 200C). This makes it possible to compare the models 
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1 Throughout this work we'll use the term double white 
dwarf instead of double degenerate, which is commonly used, 



to the observations in more detail. 

In this paper we present a new population synthesis 
for double white dwarfs, which is different from previous 
studies in three aspects. The first are some differences in 
the modelling of the binary evolution, in particular the 
description of a common envelope without spiral-in, in 
which the change in orbit is governed by conservation of 
angular momentum, rather than of energy (Sect. ^). The 
second new aspect is the use of detailed models for the 



because the term double degenerate is sometimes used for 
white dwarf - neutron star or double neutron star binaries. 
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cooling of white dwarfs (Sect. 4.3), which are important 
because it is the rate of cooling which to a large extent 
determines how long a white dwarf remains detectable in 
a magnitude-limited observed sample. The third new as- 
pect is that we use different models of the star forma- 
tion history (Sect. 0). Results are presented in Sect. [|and 
discussed in Sect. M The conclusions are summarised in 
Sect. ^. In the Appendix some details of our population 
synthesis are described . 



2. Binary and single star evolution; the formation 
of double white dwarfs 



The code we use is based on the code described by Porte- 



gies Zwart & Vcrbunt (1996) and Portegies Zwart & Yun 



gelson (1998), but has been modified in two respects; the 
white dwarf masses and the treatment of unstable mass 
transfer. 

2.1. White dwarf masses 

The masses of white dwarfs in binaries provide important 
observational constraints on evolution models. Therefore 
we have improved the treatment of the formation of white 
dwarfs in our binary evolution models by keeping more ac- 
curate track of the gro wth of the mass of the core. Details 
are given in Appendix |A.l.l| . 

2.2. Unstable mass transfer 

There exist two "standard" scenarios for the formation of 
close double white dwarfs. In the first, the binary experi- 
ences two stages of unstable mass transfer in which a com- 
mon envelope is formed. The change of the binary orbital 
separation in a common envelope is treated on the base of 
a balance between orbital energy and the binding energy 



of the envelope of the mass-losing star ( 


Paczyhski 1976; 


Tutukov & Yungelson 1979; 


Webbink 1984; 


Ibcn & Livio 



1993). The second scenario assumes that the first-born 



white dwarf of the pair is formed via stable mass transfer, 
like in Algol-type binaries (possibly accompanied by some 
loss of mass and angular momentum from the system) and 
the second white dwarf is formed via a common envelope. 

Reconstruction of the evolution of three double helium 
white dwarfs with known masses of both components led 
us to the conclusion that a spiral-in could be avoided in 



the first phase of unstable mass transfer (Nelemans et al 



200C). Briefly, when the mass ratio of two stars entering 
a common envelope is not too far from unity, we assume 
that the envelope of the evolving giant is ejected without a 
spiral-in, and that the change in orbital separation is gov- 
erned by conservation of angular momentum (the equation 



used is given in Appendix A. 2. 3). We parametrise the loss 



of angular momentum from the binary with a factor 7. 
If the mass ratio is more extreme, the common envelope 
leads to a spiral-in, which is governed by the conservation 



of energy (the equation used is given in Appendix A. 2. 2). 
The efficiency with which the energy of the binary orbit 
is used to expell the envelope of the giant is parametrised 
by a factor a co A. We switch between the two descriptions 
at the mass ratio where both give the same change of the 



separation (roughly at 0.2). Nelemans et al. (2000) find 
that values of 7 = 1.75 and a CG A = 2 give the best agree- 
ment of evolution models with the observed parameters of 
three binaries in which the masses of both white dwarfs 
are known, and therefore we use these values in our cal- 
culations. 

Another novelty is what we suggest to call "double 
spiral-in". It describes the situation when the primary 
fills its Roche lobe at the time that its companion has 
also evolved off the main sequence. This kind of evolution 
can only take place when the initial mass ratio is close 
to unity. Such a mass transfer phase has hitherto been 
described with the standard common envelope formalism; 
in the same way as when the companion is still a main 
sequence star. However, if the companion is evolved, one 
might as well argue that the envelope of the smaller star 
becomes part of the common envelope, and the envelopes 
of both stars will be expelled. We propose to use the energy 
balance here, since the double core binary will in general 
not have enough angular momentum to force the envelope 
into co-rotation. An equation for the change in orbital 
separation in the case of a "double spiral-in" is derived in 
Appendix A. 2. 4 exactly analogous to the usual common 
envelope formalism (e.g. Webbink 1984 ). 



2.3. Examples 

Before discussing effects that influence the double white 
dwarf population as a whole we discuss some typical ex- 
amples of binary evolution leading to close double white 
dwarfs, to illustrate some of the assumptions used in our 
models. For details of the treatment of binary evolution 



Verbunt (1996) and the 



we refer to Portegies Zwart 
Appendix. 



2.3.1. Double helium white dwarfs 

The most common double white dwarfs consist of two he- 
lium white dwarfs (Sect. |6.l| ). These white dwarfs descend 
from systems in which both stars have M 5s 2.3 Mq and 
fill their Roche lobes before He ignition in their degen- 
erate cores. In Fig. [j] (top left) we show an example of 
the formation of such a system. We start with a binary 
with an orbital period of 40 days and components of 1.4 
and 1.1 Mq. The primary fills its Roche lobe after 3 Gyrs, 
at which moment it has already evolved up the first giant 
branch and has lost ^0.13 Mq in a stellar wind. When the 
star fills its Roche lobe it has a deep convective envelope, 
so the mass transfer is unstable. We apply the envelope 
ejection formalism to d escrib e the mass transfer with a 
7- value of 1.75 (see Eq. A. 16] ). The core of the donor be- 
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Fig. 1. Evolutionary scenarios for the formation of a double helium white dwarf (top left), a double CO white dwarf 
(top right) and the CO+He and He+CO pairs (bottom ones). Note that t he sc ales in the panels differ as indicated by 
the 100 Rq rulers at the bottom. For a more detailed discussion see Sect. 2.3 



comes a 0.31 Mq helium white dwarf. The orbital period rather extreme and a common envelope is formed in which 
of the system hardly changes. After 4 Gyr, when the first the orbit shrinks dramatically, 
formed white dwarf has already cooled to very low lumi- 
nosity, the secondary fills its Roche lobe and has a deep 
convective envelope. Mass loss again proceeds on dynam- 
ical time scale, but the mass ratio of the components is 
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2.3.2. Double CO white dwarfs 

Most double CO white dwarfs are formed in systems which 
are initially so wide that both mass transfer phases take 
place when the star is on the AGB and its core consists 
already of CO, such that a CO white dwarfs are formed 
directly. An example is shown in Fig. [l] (top right). In 
the first phase of mass transfer the change of the orbital 
separation is regulated by the conservation of angular mo- 
mentum during envelope ejection, according to equation 
(Eq. A. 16), while in the second phase of mass transfer 



Table 1. Models and their parameters. IMF is always 
according to |Miller fc Scalo (1979| ). The SFR is either ex- 
ponentially decaying (Eq. [4|) or constant. The column "% 
binaries" gives the initial binary fraction in the popula- 
tion, the column "cooling" gives the cooling model (see 
Sect. IO) 



spiral- in is described by Eq. (A. 14). 

Much less frequently, CO white dwarfs are formed by 
stars more massive than 2.3 Mq which fill their Roche 
lobe when they have a nondegenerate core, before helium 
ignition. Roche lobe overflow then results in the forma- 
tion of a low-mass helium star. A brief additional phase of 
mass transfer may happen, if the helium star expands to 
giant dimensions during helium shell burning. This is the 
case for 0.8 <; Af H c/M <3 (see Appendix |A.1.2j ). After 
exhaustion of helium in its core, the helium star becomes 
a CO white dwarf. 

2.3.3. CO white dwarfs with He companions 

In Fig. (bottom left) we show an example in which the 
CO white dwarf is formed first. It starts with a more ex- 
treme mass ratio and a relatively wide orbit, which shrinks 
in a phase of envelope ejection. The secondary does not 
accrete anything and fills its Roche lobe when it ascends 
the first giant branch, having a degenerate helium core. It 
then evolves into a helium white dwarf. 

In t he oocond example (ahown in Fig. [l|; bottom right), 
the sy atom ovolvco through a otablo maoo exchange phaoo 
because the primary has a radiative envelope when it fills 
its Roche lobe. Part of the transferred mass is lost from 
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initialise 250,000 "zero-age" binaries and evolve these bi- 
naries according to simplified prescriptions for single and 
binary star evolution, including stellar wind, mass transfer 
(which may involve loss of mass and angular momentum 
from the binary), common envelopes and supernovae. 

For each initial binary the mass Af; of the more massive 
component, the mass ratio = m\jM\ < 1, where mi 
is the mass of the less massive component, the orbital 
separation a; and eccentricity e\ are chosen randomly from 
distributions given by 

Prob(Afi) MS79 for 0.96 Mq < Mi < 11 Mq, 
Prob(<7;) oc const, for < (ft < 1, 



Prob(a;) cx 



for < logai/i? < 6, 



(1) 



Prob(ei) cx 2e; for < e A < 1. 



For the primary mass we use the approximation of Egglc- 



the systems (see Appendix A.2.1 ). The orbit widens and 
the primary forms a helium white dwarf when it has trans- 
ferred all its envelope to its companion. The secondary ac- 
cretes so much mass that it becomes too massive to form a 
helium white dwarf. The secondary fills its Roche lobe on 
the AGB to form a CO white dwarf in a common envelope 
in which the orbital separation reduces strongly. Because 
of the differential cooling (Sec. 4.3) the CO white dwarf, 
despite the fact that it is formed last, can become fainter 
than its helium companion. Since the probability to fill 
their Roche lobe when the star has a radiative envelope, 
is low for low-mass stars, the scenario in which the helium 
white dwarf is formed first is less likely (see Sect. |6|). 

3. A model for the current population of white 
dwarfs in the Galaxy 

We model the current population of double and single 
white dwarfs in the Galaxy using population synthesis and 
compare our models with the observed populations. We 



ton et al. (1989j ) to the [Miller fc Scalo (1979| ) IMF indi- 
cated as MS79. A primary at the lower mass limit has a 
main sequence life time equal to our choice of the age of 
the Galactic disk (10 Gyr). The lower mass of less mas- 
sive component is set to 0.08 Mq, the minimum mass for 
hydrogen core burning. The distribution over separation 
is truncated at the lower end by the separation at which 
the ZAMS binary would be semi-detached. 

To investigate the effects of different cooling models 
(Sect. |4.3| ) and different assumptions about the star for- 
mation history (Sect. ^) different models have been com- 
puted (Table jl|). 

4. Modelling the observable population; white 
dwarf cooling 

To model the observable population we have to take or- 
bital evolution and selection effects into account. 



4-1- Orbital evolution of double white dwarfs 

The most important effect of orbital evolution, which is 
taken into account also in all previous studies of close bi- 
nary white dwarfs, is the disappearance from the sample 
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of the tightest systems as they merge, due to the loss of 
angular momentum via gravitational wave radiation. For 
example an 0.6 Mq+ 0.6 Mq white dwarf pair with orbital 
period of 1 hour merges in 3 x 10 7 yr. If it is located at a 



velopes is terminated by Roche lobe overflow. The amount 
of hydrogen that is left on the white dwarf depends on the 
details of this process. Fully fledged evolutionary calcula- 



distancc of 100 pc from the Sun it will disappear abruptly 
from a magnitude limited sample by merging |^| before the 
white dwarfs have become undetectable due to cooling. 



4-2. Selection effects 

The observed double white dwarfs are a biased sample. 
First, they were mainly selected for study because of their 
supp osed low mass, since this is a clear ind ication of bina- 
rity ( |Saffer et al. 198$ [Marsh et al. 1995| ). Secondly, for 
the mass determinations and the measurement of the ra- 
dial velocities the white dwarfs must be sufficiently bright. 
A third requirement is that the radial velocities must be 
large enough that they can be found, but small enough 
that s pectral lines don't get smeared out during the in 



tions of the formation of helium white dwarfs, e.g. Gian- 
none fe Giannuzzi (1970|); Sarna et al. (2000), as well as 



calculations that mimic Roche lobe overflow by mass loss 
at fixed constant rate (Driebe et al. 1998), find that the 
thickness of the residual envelope around the white dwarf 
is increasing with decreasing white dwarf mass. As a result 
the brightness at fixed age decreases mon oton ically with 
increasing white dwarf mass (see also Fig. |A.2| ). 

However, it is not clear that these calculations are valid 
for white dwarfs formed in a common envelope. In addi- 
tion, white dwarfs may lose mass by stellar wind when 
they still have a high luminosity. Such winds are observed 
for nuclei of planetary nebulae and post-novae and could 
also be expected for He white dwarfs. Finally, white dwarfs 
with masses between ^0.2 and ~0.3 Mq experience ther- 
mal flashes (|Kippcnhahii et al. 1968 ; Webbink 1975 ; [ben 



tegratibn. |Maxtcd fc Marsh (1999|) discuss this last re- fc Tutukov 1986b|; priebe et al. 1999|; |Sarna et al. 2000| ), in 



quircmcnt in detail. Following them, we include a detec- 
tion probability in the model assuming that double white 
dwarfs in the orbital period range between 0.15 hr and 
8 . 5 da y will be detected with 100% probability and that 



which the envelopes expand. This may lead to additional 
mass loss in a temporary common envelope, especially in 
the closest systems with separations < IRq. M ass loss 
may result in extinguishing of hydrogen burning (Iben & 



above 6-5 day the detection probability decreases linearly Tutukov 1986b|; |Sarna et al. 2000[ ) 



from 1 at 8.5 d ays to at ^35 days (see Fig. 1 in Maxted 
fc Marsh 1991 ) . 



The second selection effect is related to the brightness 
of the white dwarfs, which is governed by their cooling 
curves. 

4-3. White dwarf cooling 



Iben & Tutukov (1985) noticed that for a 0.6 M white 
dwarf the maximum probability of discovery corresponds 
to a cooling age of ~ 10 8 yr. In absence of detailed cooling 
curves for low-mass white dwarfs, it was hitherto assumed 
in population synthesis studies that white dwarfs remain 
bright enough to be observed during 10 s yr, irrespective of 



Hansen (tyyy) argues that the details of the loss of the 
hydrogen envelope are very uncertain and assumes that all 
white dwarfs have a hydrogen envelope of the same mass. 
He finds that helium white dwarfs cool slower than the CO 
white dwarfs, but inside these groups, the more massive 
white dwarfs cool the slowest. The difference within the 
groups are small. 

We conclude that the cooling models are still quite 
uncertain, so we will investigate the result of assuming 
different cooling models in our population synthesis. 

The first model we compute (Al; see Table |Tj for a 
list of all computed models) uses the cooling curves as 
given by |Blocker (1995| ) for CO white dwarfs an d |DSBH9S 
for He white dwarfs as detailed in Appendix A. 1.5. For 



their mass. However, recent computations (Blocker 1995; 
Driebe et al. |1998| , hereafter pSBH98| ; Hansen |1999| ) indi- 
cate that helium white dwarfs cool more slowly than CO 
white dwarfs, for two reasons. First, helium cores contain 
a higher number of ions than carbon-oxygen cores of the 
same mas s, they store more heat and are brighter at the 
same age ( Hansen 1999 ). Second, if the mass of the hydro- 
gen envelope of the white dwarf exceeds a critical value, 
pp-reactions remain the main source of energy do wn to 
effective temperatures well below 10 4 K (Webbink 1975 ; 
DSBH98| ; Sarna et al. |2000| ). This residual burning may 
lead to a significant slow-down of the cooling. 

White dwarfs in close binaries form when the evolution 
of (sub)giants with degenerate cores and hydrogen-rich en- 



2 Note, however, that just before mergi ng: white dwarfs m ay 
become quite bright due to tidal heating (Iben et al. 1998). 



the second model (A2) we made a crude estimate of the 
cooling curves for the case that the thermal flashes or a 
stellar wind reduce the mass of the hydrogen envelope and 
terminate the residual burning of hydrogen. We apply this 
to white dwarfs with masses below 0.3 Mq, and model 
all these white dwarfs identically and simply with cooling 
curves for a more massive (faster cooling) white dwarf of 
0.46 Mq. To compare with the previous investigators, we 
include one model (A3) in which all white dwarfs can be 
se en for 100 Myr s. We did not model the cooling curves 
of Hansen (1999] ), because no data for L > O.OILq are 
given. 



4-4- Magnitude limited samples and local space densities 

To convert the total Galactic population to a local pop- 
ulation and to compute a magnitude limited sample, we 
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assume a distribution of all single and binary stars in the 
galactic disk of the form 



p(R, z) = po e~ R l H sech(z/h) 2 



pc 



(2) 



where we use H — 2.5 kpc ( fjackett 1997 ) and h = 200 pc, 
neglecting the age and mass dependence of h. 

To construct a magnitude limited sample, we com- 
pute the magnitude for all model systems from the cooling 
curves and estimate the contribution of each model sys- 
tem from Eq. The absolute visual magnitudes along 
the cooling curves are derived using bolometric corrections 
after |Eggleton et al. (1989[ ). 

From Eq. (|2j) the local (R = 8.5kpc, z = 30pc) space 
density {p i q) of any type of system is related to the total 
number in the Galaxy (Nj) by: 



p itQ = A,/4.8 x 10 1 



pc 



(3) 



5. Star formation history 



Some progenitors of white dwarfs are formed long ago. 
Therefore the history of star formation in the Galaxy af- 
fects the contribution of old stars to the population of local 
white dwarfs. To study this we compute different models. 

For models A and B (see Table [l]), we model the star 
formation history of the galactic disk as 

SFR(i) = 15 exp(-*/r) M yr" 1 (4) 

where t = 7 Gyr. It gives a current rate of 3.6MQyr" 



which is compatible with observational estimates (Rana 



199l| ; |van den Hock fe de Jong 1997| ). The integrated SFR 



i.e. the amount of matter that has been turned into stars 
over the whole history of the galactic disk (10 Gyr) with 
this equation is ^8 x 10 10 Mq which is higher than the 
current mass of the disk, since part of the gas that is 
turned into stars is given back to the ISM by supernovae 
and stellar winds. 

For models C and D we use a constant SFR of 4Mq 



yr (as Tutukov & Yungelson 1993). We use an age of 



the disk of 10 Gyr, while Tutukov fc Yungelson (1993| ) use 
15 Gyr. Model D also allows us to compare our results 
with previous studies ( |ITY97| and |HAN98| ; see Sect. |). 

Most binary population synthesis calculations take a 
binary fraction of 100%. Since we want to compare our 
models with the observed fraction of close double white 
dwarfs among all white dwarfs, we present models with 
100% binaries (models B and D); and with 50% binaries 
and 50% single stars, i.e. with 2/3 of all stars in binaries 
(models A and C). 

6. Results 

Our results are presented in the next five subsections. In 



Sect. 6.1 we give the birth rates and total number of dou- 



ble white dwarfs in the Galaxy. These numbers allow a 
detailed comparison with results of earlier studies, which 



we defer to Sect. 0. They cannot be compared with ob- 
servations directly, with the exception of the SN la rate. 
For comparison with the observed sample, described in 



Sect. S.2, we compute magnitude limited samples in the 



remaining sections. In Sect. 3.3 the distribution over peri 



ods and masses is compared with the observations, which 
constrains the cooling models. Comparison of the mass 
ratio distribution with the observations gives further sup- 
port for our new d escr iption of a co mmon envelope with- 
out spiral-in (Sect. |6.4[) . In Sect. 6J) we compare our model 
with the tota l po pulation of single and binary white dwarfs 
and in Sect. |6.6| we compare models that differ in the as- 
sumed star formation history with the observed rate of 
PN formation and the local space density of white dwarfs. 

6.1. Birth rates and numbers 

In Table || the birth rates for all models are given. Ac- 
cording to Eq. (Q) the mass of a binary is on average 1.5 
times the mass of a single star. For each binary in models 
A and C we also form a single star, i.e. per binary a total 
of 2.5 times the mass of a single star is formed (1.5 for the 
binary, 1 for the single star). For models B and D only 
1.5 times the mass of a single star is formed per binary. 
Thus for the same SFR in MQyr -1 the frequency of each 
process involving a binary of the models A and C is 0.6 
times that in models B and D. 

For model A the current birth rate for close double 
white dwarfs is 4.8 x 10~ 2 yr -1 in the Galaxy. The ex- 
pected total population of close binary white dwarfs in the 
galactic disk is ~ 2.5 x 10 s (see Table ||). 

The double white dwarfs are of the following types: 
53% contains two helium white dwarfs; 25% two CO white 
dwarfs; in 14% a CO white dwarf is formed first and a 
helium white dwarf later and in 6% a helium white dwarf 
is formed followed by the formation of a CO white dwarf. 
The remaining 1% of the double white dwarfs contains 
an ONeMg white dwarf. The CO white dwarfs can be so 
called hybrid white dwarfs; having CO cores and thick 



helium envelopes ( |Iben fc Tutukov 1985| , [19871 ). of the 
double CO white dwarfs, 6% contains one and 5% two 
hybrid white dwarfs. In the mixed pairs the CO white 
dwarf is a hybrid in 20% of the cases. 

Forty eight percent of all systems are close enough to 
be brought into contact within a Hubble time. Most are 
expected to merge. The estimated current merger rate of 
white dwarfs is 2.2 x 10 -2 yr -1 . The current merger rate 
of pairs that have a total mass larger than the Chan- 
drasekhar limit (M ch = 1.44M Q ) is 3.2 x lO^yr" 1 . 
Since the merging of binary CO white dwarfs with a 
combined mass in excess of Mch is a viable model for 



type la SNe (sec Livio 1999, for the most recent review), 



our model rate can be compared with the SN la rate of 
~(4±1) x 10~ 3 yr~ 1 for Sbc type galaxies like our own 
( Cappellaro et al. 1999 ). In 19% of the systems that come 
into contact the ensuing mass transfer is stable and an 
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Model 


SFH 


% bin 


"(wd,wd) 


"merge 


SN la 


^AMCVn 


#(wd, wd) 








(io- 2 ) 


(IO" 2 ) 


(IO" 3 ) 


(IO" 3 ) 


(10 8 ) 


A 


Exp 


50 


4.8 


2.2 


3.2 


4.6 


2.5 


B 


Exp 


100 


8.1 


3.6 


5.4 


7.8 


4.1 


C 


Cnst 


50 


3.2 


1.6 


3.4 


3.1 


1.2 


D 


Cnst 


100 


5.3 


2.8 


5.8 


5.2 


1.9 


[TY97 1 


Cnst 


100 


8.7 


2.4 


2.7 


12.0 


3.5 


HAN98 


1 Cnst 


100 


3.2 


3.1 


2.9 


26 


1.0 



1 Note that [TY97 and HAN98 used a normalisation that is higher than we 
use for model D by factors ~1.4 and ~1.1 respectively (see Sect. [r.l| ). 



Table 2. Birth and event rates and 
numbers for the different models. All 



of yr 1 in the Galaxy. All numbers 
(0) are total numbers in the Galaxy. 



Close double white dwarfs are repre- 
sented with (wd, wd). See Sect. 6.1 for 
a discussion of these rates. For com- 
parison (in Sect. [7.1|) we also include 



■numbers computed by the code from 



[TY97| but using an age of the galactic 



disk of 10 Gyr instead of the 15 Gyr 
used by [TY97; and numbers of model 



1 of HAN98 



Table 3. Parameters of known close double white dwarfs 
(first 14 entries) and subdwarfs with white dwarf com- 
panions, m denotes the mass of the visible white dwarf or 
subdwarf. The mass ratio q is defined as the mass of the 
brighter star of the pair over the mass of the companion 
For ref erences see [Maxtcd k Marsh (1999[); [Moran ct al 



(1999|)] |Marsh (1999j ); and |Maxted et al. (200C|). The mass 



of 0136+768 is corrected for a misprint in |Maxtcd k Marsh 



(1999) 
(1997]) 
from 



for 0135+052 the new mass given in [Bergeron ct al 
is taken 



Data for th e sdB star KPD 0422+5421 are 
Orosz k Wade (1999| ) and for KPD 1930+2752 from 



Maxted et al. (2000). The remaining sdB stars do not have 
reliable mass estimates. 



WD/sdB 


P(d) 


q 


m 


sdB 


P(d) 


0135-052 


1.556 


0.90 


0.25 


0101+039 


0.570 


0136+768 


1.407 


1.31 


0.44 


0940+068 


8.33 


0957-666 


0.061 


1.14 


0.37 


1101+249 


0.354 


1022+050 


1.157 




0.35 


1432+159 


0.225 


1101+364 


0.145 


0.87 


0.31 


1538+269 


2.50 


1202+608 


1.493 




0.40 


2345+318 


0.241 


1204+450 


1.603 


1.00 


0.51 






1241-010 


3.347 




0.31 






1317+453 


4.872 




0.33 






1704+481A 


0.145 


0.7 


0.39 






1713+332 


1.123 




0.38 






1824+040 


6.266 




0.39 






2032+188 


5.084 




0.36 






2331+290 


0.167 




0.39 






KPD 0422+5421 


0.090 


0.96 


0.51 






KPD 1930+2752 


0.095 


0.52 


0.5 







6.2. Observed sample of double white dwarfs 

The properties of the observed double white dwarfs with 
which we will compare our models are summarised in Ta- 
ble |. Only WD 1204+450 and WD 1704+481 are likely to 
contain CO white dwarfs, having components with masses 
higher than 0.46 M^; the limitin g mass to form a helium 
white dwarf (Sweigart et al. 199C). The remaining systems 
are probably helium white dwarfs. In principle in the mass 
range M ~ 0.35 — 0.45 Mq white dwarfs could also be hy- 
brid; However in this range the probability for a white 
dwarf to be hybrid is 4 - 5 times lower than to be a he- 
lium white dwarf, because hybrid white dwarfs originate 
from more massive stars which fill their Roche lobe in a 
narrow period range (see, ho wever, an example of s uch a 
scenario for WD 0957-666 in |Nclcmans ct al. 2000|) . We 
assume 0.05 Mq for the uncertainty in the estimates of 
the masses of white dwarfs, which may be somewhat op- 
timistic. 

Table ^ also includes data on subdwarf B stars with 
suspected white dwarf companions. Subdwarf B (sdB) 
stars are hot, helium rich objects which are thought to 
be helium burning remnants of stars which lost their hy- 
drogen envelope. When their helium burning has stopped 
they will become white dwarfs. Of special interest are 
KPD 0422+5421 ( |Koen et al. 1998j; |Orosz k Wade 1999Q 
and KPD1930+2752 ([Maxted et al. 2000|). With orbital 



periods as short as 0.09 and 0.095 days, respectively, their 
components will inevitably merge. In both systems the 
sdB components will become white dwarfs before the stars 
merge. In KPD1930+2752 the total mass of the compo- 
nents is close to the Chandrasckhar mass or even exceeds 
it. That makes this system the only currently known can- 
didate progenitor for a SN la. 



interacting double white dwarf (identified with AM CVn 
stars) is formed. The model birth rate of AM CVn systems 
is 4.6 x 10" 3 yr" 1 (see Table f|. 



6.3. Period - mass distribution; constraints on cooling 
models 

The observed quantities that are determined for all double 
white dwarfs are the orbital period and the mass of the 
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Fig. 2. Model population of double white dwarfs as function of orbital period and mass of the brighter white dwarf of 
the pair. Top left: distribution of the double white dwarfs that are currently born for models A. This is independent 
of cooling. In the remaining three plots we show t he current ly vi sible populatio n of double white dwarfs for different 
cooling mod els: (top right) cooling according to DSBH98 and Blocker (1995 , model Al); (bottom right) cooling 
according to DSBH98 , but with faster cooling for WD with masses below 0.3 Mq (model A2). Both plots are for a 
limiting magnitude V\i m = 15; (bottom left) with constant cooling time of 100 Myr (model A3, note that in this case 
we only obtain the total number of potentially visible double white dwarfs in the Galaxy and we can not construct a 
magnitude limited sample). For comparison, we also plot the observed binary white dwarfs. 



brighter white dwarf. Following gaffer et al. (199lf ), we 
plot in Fig. H the P or b — m distributions of the frequency 
of occurrence for the white dwarfs which are born at this 
moment and for the simulated magnitude limited sample 
for the models with different cooling prescriptions, (mod- 
els Al, A2 and A3; see Table 0), where we assume Vii m 
= 15 as the limiting magnitude of the sample^. For m we 
always use the mass of the brighter white dwarf. In gen- 
eral the brighter white dwarf is the one that was formed 
last, but occasionally, it is the one that was formed first 
as explained in Sect. 2.3.3 . For comparison, we also plot 
the observed binary white dwarfs in Fig. ||. 

There is a clear correlation between the mass of new- 
born low-mass (He) white dwarf and the orbital period of 
the pair. This can be understood as a consequence of the 



3 The P — m distribution does not qualitatively change if we 
increase V\i m by one or two magnitudes, since we still deal with 
very nearby objects. 



existence of a steep core mass - radius relation for giant s 
with degenerate helium cores ( Refsdal fc Weigert 1970 ). 
Giants with more massive cores (forming more massive 
white dwarfs) have much larger radii and thus smaller 
binding energies. To expell the envelope in the common 
envelope, less orbital energy has to be used, leading to 
a larger orbital period. The spread in the distribution is 
caused by the difference in the masses of the progenitors 
and different companion masses. 

In the simulated population of binary white dwarfs 
there are three distinct groups of stars: He dwarfs with 
masses below 0.45 Mq, hybrid white dwarfs with masses 
in majority between 0.4 and 0.5 Mq and periods around 
a few hours, and CO ones with masses above 0.5 Mq. The 
last groups are clearly dominated by the lowest mass ob- 
jects. The lowest mass CO white dwarfs are descendants of 
most numerous initial binaries with masses of components 
1 - 2M . 
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Fig. 3. Cumulative dis tribution of periods. Solid line for 
our best model (A2); DSBH9S cooling, but with lower 



luminosity due to thermal flashes for white dwa rfs with 
masses below 0.3 Mq. Dashed line for IDSBH98 without 



modifications (model Al) and dash dotted line for con- 
stant cooling time of 100 Myr (model A3). Open squares 
for the observed double white dwarfs, filled circles give the 
observed systems including the sdB binaries (Table ||). 

The different cooling models result in very different 
predicted observable distributions. Model Al where the 



cooling curves of DSBH98 are applied favours low mass 
white dwarfs to such an extent that almost all observed 
white dwarfs are expected to have masses below 0.3 Mq. 
This is in clear contrast with the observations, in which 
all but one white dwarf have a mass above 0.3 Mq. Re- 
duced cooling times for white dwarfs with masses below 
0.3 Mq (model A2) improves this situation. Model A3, 
with a constant cooling time (so essentially only affected 
by merging due to GWR), seems to fit all observed sys- 
tems also nicely. However, a complementary comparison 
with the observations as given by cumulative distributions 
of the periods (Fig. ||), shows that model A2 fits the data 
best, and that model A3 predicts too many short period 
systems. 

The observed period distribution for double white 
dwarfs shows a gap between 0.5 and 1 day, which is not 
present in our models. If we include also sdB binaries, the 
gap is partially filled in. More systems must be found to 
determine whether the gap is real. 

The comparison of our models with observations sug- 
gests that white dwarfs with masses below 0.3 Mq cool 



faster than predicted by DSBH98. Mass loss in thermal 



flashes and a stellar wind may be the cause of this. 

The model sample of detectable systems is totally dom- 
inated by He white dwarfs with long cooling times. Given 
our model birth rates and the cooling curves we apply, 
we estimate the number of double white dwarfs to be de- 
tected in a sample limited by Vum = 15 as 220 of which 
only 10 are CO white dwarfs for model A2. Roughly one of 



//'/ P 
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-2-10 1 
log P (days) 

Fig. 4. Cumulative distribution of periods. Solid line for 
model A2 as in Fig. dashed line for the same model but 
with a cc X = 1, dash-dotted line for a model with 7 = 1.5 
and finally the dotted line for model C (constant SFR). 

Table 4. Number of observable white dwarfs, close double 
white dwarfs and SN la progenitors as function of the 
limiting magnitude of the sample for model A2. 





#wd 


#wdwd 


#SN la prog 


15.0 


855 


220 


0.9 


15.5 


1789 


421 


1.7 


16.0 


3661 


789 


3.2 


17.0 


12155 


2551 


11.2 



these is expected to merge within a Hubble time having a 
total mass above Mch- For future observations we give in 
Table |^ a list of expected number of systems for different 
limiting magnitudes. 

It should be noted that these numbers are uncertain. 
This is illustrated by the range in birth rates for the differ- 
ent models (Table 0) and by the differences with previous 
studies (see Sect. |7.1| ). Additional uncertainties are intro- 
duced by our limited knowledge of the initial distributions 
(Eq. |l|) and the uncertainties in th e cooling and the Galac- 
tic model (Eq. (||)). For example Yungelson et al. (1994| ) 
compare models with two different qi distributions (one 
peaked towards gi ~ 1) and show that the birth rates dif- 
fer by a factor ~ 1.7. In general the relative statistics of 
the model are more reliable than the absolute statistics. 

Before turning to the mass ratio distribution, we il- 
lustrate the influence of the model parameters we choose. 
We do this by showing cumulative period distributions for 
some models with different parameters in Fig. H a cc A = 
1 (dashed line) and 7 = 1.5 (dash-dotted line). It shows 
that the change in parameters influences the distributions 
less than the different cooling models discussed above, al- 
though the observations favour a higher a cc A. We also 
included the cumulative distribution for model C (with 
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Fig. 5. Top: current population of double white dwarfs 
as function of orbital period and mass ratio, for model 
A2, a limiting magnitude of 15 and a maximal ratio of 
luminosities of 5. Bottom: the same for a run in which 
the first phase of mass transfer is treat ed a s a stan dard 
common envelope, as is done by 1TY97 and HAN98 . For 
comparison, we also plot the observed binary white dwarfs. 



a constant SFR; dotted line) which differs from that for 
model A2 in that it has fewer long period systems. This 
is a consequence of the larger relative importance of old, 
low-mass progenitor binaries in model A2, which lose less 
mass and thus shrink less in the first phase of mass trans- 
fer (see Eq. A. 16 ). 



6.4- Period - mass ratio distribution 

Our assumption that a common envelope can be avoided 
in the first phase of mass transfer between a giant and 
a main-sequence star, is reflected in the mass ratios of 
the model systems. A clear prediction of the model is 
that close binary white dwarfs must concentrate to q = 
m/M ~ 1. For the observed systems, the mass ratio can 
only be determined if both components can be seen which 
in practice requires that the luminosity of the fainter com- 
ponent is more than 20% of that of the brighter component 




2.5 



m/M 



Fig. 6. Cumulative mass ratio distributions for the mod- 
els A 2 (s olid line) and A' (dotted line) as explained in 
Sect. 14. The observed mass ratio's are plotted as the 
open squares. 



(Moran et al. 200C). Applying this selection criterium to 
the theoretical model, we obtain the distribution shown 
m Fig. | for the magnitude limited sample. Note that 
since lower mass white dwarfs cool slower this selection 
criterium favours systems with mass ratios above unity. 
In the same figure we also show the observed systems. 

For comparison we also computed a run (A') in which 
we used the standard common envelope treatment for the 



first phase of mass transfer, which is done by [TY97 and 



HAN98. The fraction of double white dwarfs for which the 



mass ratio can be determined according to the selection 
criterium of a luminosity ratio greater than 0.2, is 27% for 
model A2 and 24% for model A'. In a total of 14 systems 
one thus expects 4±2 and 3±2 systems of which the mass 
ratio can be determined. Model A2 fits this number better, 
but the numbers are too small to draw conclusions. The 
distribution of mass ratios in model A' (Fig. |5|, bottom) 
however clearly does not describe the observations as well 
as our model A2, as illustrated in more detail in a plot 
where the cumulative mass ratio distributions of the two 
models and the observations are shown (Fig. ||). 

6. 5. Mass spectrum of the white dwarf population; 
constraints on the binary fraction 

Figs. |l] and || show the model spectrum of white dwarf 
masses for models B and A2, including both single and 
double white dwarfs for a limiting magnitude V\i m = 15. 
For this plot we consider as "single" white dwarfs all ob- 
jects that were born in initially wide pairs, single merger 
products, white dwarfs that became single as a result of 
binary disruption by SN explosions, white dwarfs in close 
pairs which are brighter than their main-sequence com- 
panions and genuine single white dwarfs for the models 
with an initial binary fraction smaller than 100%. 
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Fig. 7. Mass spectrum of all white dwarfs for model 
B (100% binaries). Members of close double white 
dwarfs are in grey. The cumulative distribution is 
shown as the solid black line. For comparison, the 
grey line shows the cumulative distribution of the ob- 
served systems (Fig. 0). 
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Fig. 9. Mass spectr um of observed whit e d warfs 
Data are taken fro m Bergeron ct al. (1992 ) and Bra- 
gaglia ;t al. (1995 ). The solid line is the cumulative 
distribution. 
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Fig. 8. Mass spectrum of all white dwarfs for model 
A2 ( initial binary fraction of 50%) Double white 
dwarfs are in grey. The cumulative distribution is 
shown as solid black line and cumulative distribution 
of observed systems as the grey line 
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Fig. 10. Mass spectrum of all white dwarfs as in 
Fig. U in a model in which all helium white dwarfs 
cool like a 0.4 Mq dwarf and all CO white dwarfs 
cool like a 0.6 Mq white dwarf. Lines are cumulative 
distributions for the model (black) and the observa- 
tions (grey) 



These model spectra can be compared to the observed 
mass spectrum of DA white dwarfs studied by Bergeron 
et al. ( [L992D and |Bragaglia et al. (1995| ), shown in Fig. |9| 
The latter distribution may have to be shifted to higher 
masses by about 0.05 Mq, if one uses models of white 
dwarfs with thick hydrogen envelopes for mass estimates 
(Napiwotzki et al. 1999). Clearly, a binary fraction of 50% 
fits the observed sample better, if indeed helium white 
dwarfs cool much slower than CO white dwarfs. We can 
also compare the absolute numbers. |Maxted & Marsh 



(1999|)] conclude that the fraction of close double white 



dwarfs among DA white dwarfs is between 1.7 and 19 % 
with 95% confidence. For model B the fraction of close 
white dwarfs is ^43% (853 white dwarfs of which 368 are 
close pairs), for model A2 is is ^26% (855 white dwarfs 
and 220 close pairs). Note that this fraction slightly de- 
creases for higher limiting magnitudes because the single 
white dwarfs are more massive and thus generally dimmer, 
sampling a different fraction of Galaxy. An even lower bi- 
nary fraction apparently would fit the data better, but is 
in conflict with the estimated fraction of binaries among 
normal main sequence binaries (Abt 1983; Duqucnnoy & 
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Table 5. Galactic number and local space density of white 
dwarfs; and Galactic and local PN formation rate for the 
models A and C. Unit of the PN formation rates is yr" 1 ; 
unit for p wdi © is pc~ 3 . The ranges of observed values are 
given for comparison. For references and discussion see 
Sect. 1611 



We can compare these numbers with the observational 
estimates for the local PN formation rate of 3 x 10 -12 



pc 3 yr 1 ( Pottasch 1996| ) and the local space density 



of white dwarfs, which range from e.g. 4.2 x 10 3 pc 3 



(Knox et al. 199£) through 7.6Tgj x 10~ 3 pc~ 3 (Oswalt 



et al. 1995) and 10 x 10~ 3 pc~ 3 (Ruiz fc Takamiya 1995) 



Model SFH % bin #wd i^>n p wd f PN q 

10 9 (1Q-' 3 ) (10- f2 ) 



A 
C 

Obs 



Exp 
Const 



50 
50 



9.2 
4.1 



1.1 
0.8 



19 

8.5 
4-20 



2.3 
1.7 

3 



Mayor 1991). However this number highly depends on un- 
certain selection effects. 

There are some features in the model mass spectrum in 
model A2 that appear to be in conflict with observations. 
The first is the clear trend that with the cooling models of 
even with our modifications, there should be an 



DSBH98 



increasing number of helium white dwarfs towards lower 
masses. The observed distribution is flat. A very simple 
numerical experiment in which we assign a cooling curve 
to all helium white dwarfs as the one for a 0.414 Mq white 
dwarf according to DSBH98j and a cooling curve a s for a 
0.605 Mq white dwarf according to Blocker (1995 ) for all 
CO white dwarfs (Fig. |l0|), shows that an equal cooling 
time for all helium white dwarfs seems to be in better 
agreement with the observations. It has a fraction of dou- 
ble white dwarfs of 18%. 

Another feature is the absence of stars with 0.45 5s 
M/Mq Ss 0.5 in the model distributions. This is a con- 
sequence of the fact that in this interval in our models 
only hybrid white dwarfs can be present, which have a 
low formation probability (see Sect. |6.2|) . 

We conclude that an initial binary fraction of 50% can 
explain the observed close binary fraction in the white 
dwarf population. The shape of the mass spectrum, espe- 
cially for the helium white dwarfs is a challenge for de- 
tailed mass determinations and cooling models. 

6. 6. Birth rate of PN and local WD space density; 
constraints on the star formation history 

Finally, we compare models A and C (see Table Q) , which 
differ only by the assumed star formation history. The 
star formation rate was probably higher in the past than 
at present and some (double) white dwarfs descend from 
stars that are formed just after the galactic disk was 
formed. 

Table | gives the formation rates of PN and the total 
number of white dwarfs in the Galaxy for models A and 
C. The total number of white dwarfs is computed by ex- 
cluding all white dwarfs in binaries where the companion 
is brighter. The local density of white dwarfs and PN rate 
are computed with Eq. (S) as described in Sect. iA. 



to 20 ± 7 x 10~ 3 pc" 3 flFcstin 1998^ 7 



This list shows the large uncertainty in the observed 
local space density of white dwarfs. It appears that the 
lower values are somewhat favoured in the literature. Both 
models A and C appear for the moment to be consistent 
with the observed local white dwarf space density and with 
the PN formation rate. However, we prefer model A2 since 
it fits the period distribution better (see Fig. Q). 

The ratio of the local space density of white dwarfs 
to the current local PN formation rate could in princi- 
ple serve as a diagnostic for the star formation history of 
the Galaxy, given better knowledge of p wdi Q, which crit- 
ically depends on the estimates of the incompleteness of 
the observed white dwarf samples and the applied cooling 
curves. 

7. Discussion: comparison with previous studies 

We now compare our work with the results of previous 



studies; in particular the most recent studies of [ben et al 
(1997| , |[TY97D and |Han (1998| , |HAN98j ). 



7.1. Birth rates 

In Table | we show the birth rates of close double white 
dwarfs for the different models. We also include numbers 



from HAN98 (model 1) and a set of numbers computed 
with the same code as used in 1TY97, but for an age of 



the galactic disk of 10 Gyr, as our models. The numbers 



of HAN98 are for an age of the disk of 15 Gyr. Our model 
D is the closest to the models of IITY97I and IHAN98L as- 



suming a constant SFR and 100% binaries. To estimate 
the influence of the binary evolution models only in com- 
paring the different models we correct for their different 
normalisations. 



In the recomputed ITY97 model the formation rate of 



interacting binaries in which the primary evolves within 
the age of the Galaxy is 0.35 yr -1 . In our model D this 
number is 0.25 yr -1 . I n the fo llowing we therefore multiply 
the formation rates of 1TY97 as given in Table ||with 0.71. 
In the model of HAN98 one binary with a primary 



mass above 0.8 Mq is formed in the Galaxy annually with 
log a; < 6.76, i.e. 0.9 binary with logai < 6, which is our 
limit to ai. Correcting for the different assumed age of the 
Galaxy we estimate this number to be 0.81; in our model 
this nu mber is .73. We thus multiply the the formation 
rates of |HAN98| as given in Table | with 0.9. 

Applying these corrections to the normalisation, we 
find that some interesting differences remain. The birth 
rate of double white dwarfs is 0.029, 0.053 and 0.062 per 
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year for |HAN98| , model D and ITY98 respectively. At the 
same time the ratio of the merger rate to the birth rate 
decreases: 0.97, 0.53 and 0.28 for these models. This can 
probably be attributed to the different treatment of the 



q ~ 0.5 — 0.7 and from HAN98] who finds typical values of 
q ~ 0.5, with a tail to q ~ 2. The difference between these 
two latter groups of models may be understood as a conse- 



commc jn envelope - HAN98 uses a common envelope spiral 



in efficiency of 1 in Webbinks (1984) formalism, while we 
use 4 (for A = 0.5, see pe Kool et al. 1987[). |ITY97j u se the 
formalism proposed by Tutukov & Yungclson (1979) with 



quence of enhanced wind in Han's model (see also Tout & 
Eggleton 1988| ), which allows wider separations before the 



an efficiency of 1. This is comparable to an efficiency of 4 - 
8 i n the We bbink formalism. This means that in the model 
of HAN98| , more systems merge in a common envelope, 
which yields a low formation rate of double white dwarfs. 
The ones that form (in general) have short periods for the 
same reason, so the ratio of merger to birth rate is high. 
In the ITY model the efficiency is higher, so more systems 
will survive both common envelopes and have generally 
wider orbits, leading to a much lower ratio of merger to 
birth rate. Our model D is somewhat in between, but also 
has the different treatment of the first mass transfer phase 
(Sect ||), in which a strong spiral-in is avoided. 

The difference between the models in the SN la rate 
(^SN ia) is related both to the total merger rate and to 
the masses of the white dwarfs. The former varies within 
a factor ~ 1.5: 0.017, 0.028, and 0.028 yr" 1 for |[TY97, 



model D, and HAN98 , while z^sn ia is higher by a factor 
2 - 3 in model D compared to the other models. This is 
caused by the initial - final mass relation in our models, 
which is derived from stellar models with core overshoot- 
ing, producing higher final masses. 

The difference in the birthrate of interacting white 
dwarfs (i^amcvh) is mainly a consequence of our treat- 
ment of the first mass transfer, which gives for mo del D 
a mass r atio dis tribu tion whi ch is peaked to 1 (Sect. 6.4), 
while in |ITY97| and |HAN98| the mass ratio is in general 
different from 1 (Sect. 7J2), favouring stable mass transfer 
and the formation of AM CVn systems. An additional fac- 
tor, which reduces the numb er of AM CVn systems is the 
assumption in model D and ITY97 that the mass trans- 
fer rate is limited by the Eddington rate. The formation 
and evolution of AM CVn stars is discussed in more detail 



Tutukov & Yungelson (1996) and Nelemans et al. (in 



preparation) 



7.2. Periods, masses and mass ratios 



Comparing our F ig. || with the corresponding Figure in 
Saffer et al. ( 1998 ), we find the same trend of higher white 



dwarf masses at longer periods. However, in our model the 
masses are higher than in the model of Saffer et al. (1998 ) 
at the same period. This is a consequence of the absence 
of a strong spiral-in in the first mass transfer phase in 
our model, contrary to the conventional common envelope 



model, as discussed in 2.2 



In our model the mass ratio distributi on is p eake d at 
q 1. This is different from the models of fTY97| and paf- 



second common envelope. The mass ratio distribution of 
our model, peaked at q ~ 1, appears to be more consistent 
with the observed mass ratio distribution. 



7.3. Cooling 

To explain the lack of observed white dwarfs with masses 
below 0.3 Mq we had to assume that these white dwarfs 



SBH98 



cool faster than predicted by the models of |T) 

The same assumption was required by [van Kerkwijk 
|st al. (2000 ), to bring the cooling age of the white dwarf 



fer et al. (1998 ) which predict a strong concentration to 



that accompanies PSR B1855+09 into agreement with the 
pulsar spin-down age; and to obtain cooling ages shorter 
than the age of the Galaxy for the white dwarfs accompa- 
nying PSR J0034-0534 and PSR J1713+0747 . 

The absence of the lowest mass white dwarfs may also 
be explained by the fact that a common envelope involving 
a giant with a low mass helium core (M c < 0.2 — Q.25M0) 
always leads to a complete merger, according to Sandquist 
|et al. (20"00| ) . However it can not explain the absence of the 
systems with 0.25 < M < 0.3Mq, which would form the 
majority of the observed systems using the full DSBH9S 
cooling (model Al; see Fig. 0). 



8. Conclusions 

We computed a model of the population of close binary 
white dwarfs and found good agreement between our 
model and the observed double white dwarf sample. A 
better agreement with observations compared to earlier 
studies is found due to two modifications. 

The first is a different treatment of unstable mass 
transfer from a giant to a main sequence star of com- 
parable mass. The second is a more detailed modelling of 
the cooling of low mass white dwarfs which became possi- 
ble because detailed evolutionary models for such white 
dwarfs became available. Our main conclusions can be 
summarised as follows. 

1. Comparing the mass distribution of the white dwarfs 
in close pairs with the observations, we find a lack of 
observed white dwarfs with masses below 0.3 Mq. This 
discrepancy can be removed with the assumption that 
low-mass white dwarfs cool faster than computed by 
Dricbe et al. (1998] ) . The same assumption removes dis- 
crepancies between observed and derived ages of low- 
mass white dwarfs that accompany recycled pulsars, 
as shown by van Kerkwijk et al. (2000[ ) . Faster cooling 
is expected if the hydrogen envelopes around low-mass 
white dwarfs are partially expelled by thermal flashes 
or a stellar wind. 
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Our models predict that the distribution of mass ra- 
tios of double white dwarfs, when corrected for obser- 



vational selection effects as described by Moran et al 
(2000D , peaks at a mass ratio of unity, consistent with 



observations. The distributions predicted in the mod- 



els by Iben et al. (1997) and Han (1998) peak at mass 
ratios of about 0.7 and above 1.5 and agree worse with 
the observations even after applying selection effects. 

3. Our models predict a distribution of orbital periods 
and masses of close double white dwarfs in satisfactory 
agreement with the observed distribution. 

4. Amongst the observed white dwarfs only a small frac- 
tion are members of a close pair. To bring our models 
into agreement with this, we have to assume an ini- 
tial binary fraction of 50% (i.e. as many single stars as 
binaries). 

5. In our models the ratio of the local number density of 
white dwarfs and the planetary nebula formation rate 
is a sensitive function of the star formation history of 
the Galaxy. Our predicted numbers are consistent with 
the observations. 

6. Using detailed cooling models we pridict that an ob- 
served sample of white dwarfs near the Sun, limited 
at the magnitude V — 15, contains 855 white dwarfs 
of which 220 are close pairs. Of these pairs only 10 
are double CO white dwarfs and only one is expected 
to merge having a combined mass above the Chan- 
drasekhar mass. The predicted merger rate in the 
Galaxy of double white dwarfs with a mass that ex- 
ceeds the Chandrasekhar mass is consistent with the 
inferre d SN la rate. 

ITY97 estimated, depending on a CG , to find one such 
pair in a sample of ~200 to ~600 white dwarfs. Revers- 
ing this argument, when the statistics become more re- 
liable, the observed number of systems with different 
types of white dwarfs could provide constraints on the 
cooling models for these white dwarfs. 
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Appendix A: Population synthesis code SeBa 

We present some changes we made to the population syn- 



thesis code SeBa (see 


Portegies Zwart & Verbunt 1996 


Portegies Zwart & Yungelson 1998 


)■ 



A.l. Stellar evolution 

As before, the treatment of stellar evolution in our code 
i s based on the fits to d etailed stellar evo lutionary models 
QEgglcton ct al. 1989| ; [Tout ct al. 1997|) , which give the 
luminosity and the radius of the stars as a function of 
time and mass. In addition to this we need the mass of 
the core and the mass loss due to stellar wind. These we 
obtain as follows. 



A. 1.1. Core mases and white dwarf masses 

For the mass of the helium core m c at the end of the main 
sequence we use (Eggleton, private communication, 1998) 

0.11 M 1 - 2 + 7 x 10~ 5 M 4 



1 + 2 x 10- 4 M 3 



(A.1) 



The mass of the core during the further evolution of the 
star is computed by integrating the growth of the core 
resulting from hydrogen shell burning: 



L 

m c = rjH — 
Jv 

where 



m = 9.6 x 1O- 12 M yr" 1 L" 1 



(A.2) 



(A.3) 



and X is the mass fraction of hydrogen in the envelope. 
During core helium burning we assume that half of the 
luminosity of the star is produced by hydrogen shell burn- 
ing, while in the double shell burning phase we assume 
that all of the luminosity is produced by the hydrogen 
shell burning. 

When giants have degenerate cores, application of a 
core mass - luminosity relation gives more accurate results 
than direct integration of the growth of the core. 

For degenerate helium cores of stars with M < 2.3Mq 
we use (Boothroyd & Sackmann 1988) 



M c = 0.146 L u 



(A.4) 



(all quantities in solar units). For degenerate CO cores of 
stars with M ^ 8Mq on the AGB we use ( proenewegci] 
|fc de Jong 1991 ) 



0.015 



L 



47488 



0.1804 



L < 15725 



M r = 0.46 



L 



M~ 



L > 15725 



(A.5) 



46818 

where the transition between the two fits occurs at M c « 
0.73Mq in stars ~ 3.5Mq where the two relations fit to- 
gether reasonably. We changed the power of the depen- 
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star and the total lifetime of the helium star. This is sug- 
gested by computations of Habets (1986| ) and gives a CO 
core of the Chandrasekhar mass for a 2.2 M helium star; 
the minimum mass to form a neutron star in our code. 

Helium stars with 0.8 ^ M ^ 3Mq expand again af- 
ter core helium exhaustion and can lose their remaining 
helium envelope in so called case BB mass transfer. The 
amount of mass that can be lost is defined as increasing 
linearly from to 45% for stars between 0.8 and 2.2 Mq 
and stays constant above 2.2 Mq. The maximum mass 
of the CO white dwarf thus formed is 1.21 Mq. Helium 
stars of lower mass (M < 0.8Mq) do not expand and re- 
tain their thick helium envelopes, forming hybrid white 



dwarfs (Iben fc Tutukov 1985) 



4.0 6.0 8.0 

ZAMS mass 



Fig. A.l. White dwarf masses as function of the ZAMS 
mass. Dashed lines are for case B mass transfer. The white 
dwarfs that descend from stars with ZAMS masses below 
2.3 Mq in case B mass transfer are helium white dwarfs. 
The two dashed lines give the minimum and maximum 
mass of the white dwarf, which depends on the orbital 
separation at the onset of the mass transfer. Solid lines are 
for case C mass transfer, which results in the formation 
of a CO white dwarf. When the ZAMS mass is above 8 
Mq the stripping of the envelope in case C mass transfer 
may prevent the formation of a neutron star, leading to 
the formation of a white dwarfs with a core consisting of 
O, Ne and Mg (shaded region). 



dence on M from -0.19 in the original paper to -0.25 be- 
cause the maximum luminosities given by our fits other- 
wise lead to white dwarf masses too high compared to 
initial - final mass relations as found from observations 



(see Groenewegen & de Jong 1993). 

The masses of CO cores formed by central He burning 
inside the helium core are defined in the same way as we 
define the relation between the mass of helium stars and 



their CO cores (see Sect. A. 1.2). 

A white dwarf forms if a component of a binary with 
M < IOMq loses its hydrogen envelope through RLOF 
either before core helium burning (case B mass transfer) 
or after helium exhaustion (case C). The masses of white 
dwarfs formed in cases B and C as function of initial mass 



are shown in Fig. A.l 



A. 1.2. Helium stars 

A helium star is formed when a star more massive than 
2.3 Mq loses its hydrogen envelope in case B mass trans- 
fer. The helium star starts core helium burning and forms 
a CO core. In our code, this core grows linearly at a rate 
given by the ratio of 65% of the initial mass of the helium 



A. 1.3. Stellar wind 

We describe mass loss in a stellar wind in a very general 
way in which the amount of wind loss increases in time 
according to 



AM W = Mi 



lost 



t + At 



'n 



(A.6) 



The exponent rj ~ 6.8 is derived from fitting stellar 
wind mass loss on the main sequence of massive stars 
(Af > 15M [Meynet et al. 1994| ), but we apply it also 
for low and intermediate mass stars. For these stars tf is 
the duration of the evolutionary phase that the star is in 



(as given by Eggleton et al. 1989). For the different evolu- 
tionary phases, the parameters Mi ost is defined as follows. 

In the Hertzsprung gap Mi ost is 1% of the total mass 
of the star. 

For the first giant branch (hydrogen she ll burning) , we 
use a fit to models of Sweigart et al. (199C) for stars with 
degenerate helium cores 



Mica* = (2.5 - M)/7.5 M 



(A.7) 



which we extend to all low and intermediate mass stars 
by setting Afi ost = above Af = 2.5M . 

On the horizontal branch Mi ost is 5% of the envelope 
mass. 

For the AGB phase we take Mi ost equal to 80% of the 
mass of the envelope of the star when it enters the early 
AGB phase. 

A. 1.4. Radii of gyration 

In the previous version of the SeBa code all gyration radii 
were set to 0.4. The gyration radius plays a role in the 
determination of the stabilit y of the mass transfer ( Porte- 
k'ics Zwart fc Vcrbunt 1996 , Appendix C.l). We now use 



the following values. 

For main-sequence st ars we use a fit to the results by 
|Clarct fc Gimcnez (199C ) . Further we classify stars either 
as radiative (stars in Hertzsprung gap and helium stars) 
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or as convective (red giants, AGB stars). A summary of 
radii of gyration are given in Table [A~l] . 

Table A.l. Gyration radii for various types of stars 
Type k 2 



Radiative stars 
Convective stars 
White dwarfs 
Neutron stars 
Black holes 



(a) punn fc Ostrikcr (1969| J 



0.03 
0.2 
0.4 
0.25" 
1 



2.0 














A. 1.5. White dwarf evolution: luminosity and radius 
We model the cooling of white dwarfs according to the 



results [of Blocker (1995 pmcf Driebe et al. J 



Luminosity 

The luminosity of white dwarfs as function of time t can 
be reasonably well modelled by 



logL 



1.41og(i/10 6 yr) 



(A. 



where L ma x is a linear fit given by 



L max = 3.83 - 4.77 M WD for 0.18 < Af WD < 0.6 (A.9) 

(mass and luminosity in solar units). Outside these limits 
L max stays constant (i.e. L max = 3 below Mwd = 0.18 and 
L max = 1 above Mwd = 0.6). For white dwarf masses 
below 0.6Mq the luminosity is constrained to be below 
logi/Lo = —0.5, for more massive white dwarfs below 



of Blocker (1995) and 



lo g L/L(V) — 2. In Fig. A. 2 we show the fits and the results 



Dricbe et al. (199£ 



Radius 

We fit ted the models of Dricbe et al. (1998 ) and Blocker 
|(1995| ), and interpolated between the fits. The fits are 
given by 



R 

Rq 



Mog(t/i0 6 yr) for M WD < O.6M . (A.10) 



The coefficients a and b are given in Table |A.2|. Fig. A. 3 



shows the fits and the corresponding detailed calculations. 
For more massive white dwarfs we use the mass-radius 



relation for zero-temperature spheres ( Nauenberg 1972 ) 



R 

Rp, 



0.01125 



'M 



WD 



Mch 



-2/3 



M 



WD 



M Ch 



2/3 



(A.ll) 



A. 1.6. Modified DSBH9S cooling 



Our modification to the cooling described above reduces 
the cooling time scale for white dwarfs with masses below 



8.0 
log t (yr) 



Fig. A. 2. White dwarf cooling tracks from Driebe et al 



1998 ) and Blocker (1995 ). Straight lines are the fits to 



these curves. The curves are for masses of 0.179, 0.300, 
0.414, 0.6 and 0.8 from top right to bottom left 




L°St (yr) 



Fig. A. 3. Wh ite dwarf radii from Driebe et al. (1998 ) and 
|B locker (1995 ). Straight lines are the fits to these curves. 
The curves are for masses of 0.179, 0.300, 0.414, 0.6 and 
0.8 from top right to bottom left 

0.3Mq. For these white dwarfs we use the cooling curve 
and the radius of a more massive, thus faster cooling white 



dwarf of 0.46 M (see Sect. 



A. 2. Mass transfer in binary stars 



As suggested by Nelemans et al. (2000] ), we distinguish 
four types of mass transfer with different outcomes: sta- 
ble mass transfer, common envelope evolution, envelope 
ejection and a double spiral-in. 
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Table A. 2. Coefficients for the fits to the white dwarf The companion does not accrete al all (see discussion in 
radii 



M WD 


a 


b 


0.2 


0.1 


0.0175 


0.4 


0.03 


0.0044 


0.6 


0.017 


0.001 


0.8 


0.011 


0.0005 



Sect. [2.2| and Nelemans et al. 2000). The change in orbital 
separation is given by 



at / Mf mf \ / Mf + rut 
~al ~ \Mi m\ ) \ Mi + m 

In this work we use 7 = 1.75. 



Mi - Mf , 
1 -7tt " • A.16) 



A. 2.1. Stable mass transfer 

The amount of mass that can be accreted by a star is 
limited by its thermal time scale 



M n 



M 



RL 
G M 



(A.12) 



If not all mass can be accreted, we assume that the ex- 
cess of mass leaves the system taking with it nj times the 
specific angular momentum of the binary. 

This assumption gives for the variation of orbital sep- 
aration 



af / Mf mf 
cti \ Mi mi 



Mf + mf 
Mi + mi 



2nj + l 



(A.13) 



We use 71 j = 2.5, which gives good agreement for the 
periods of low-m ass Algols and Be X-ray binaries ( Porte 
gies Zwart 1996|) 



A. 2. 2. Standard common envelope 

When the mass transfer is unstable due to a tidal in- 
stability, the accretor is a compact object, or the en- 
velope ejection equation gives a smaller orbital separa- 
tion, we apply the standard common envelope equation 
AE olh QWebbink 1983): 



Mi (Mi - M f ) 



Mf TO Mi TO 

2 af 2 Oi 



(A.14) 



where a cc is an efficiency parameter and A a parameter 
describing the strucure of the envelope of the giant. Both 
are uncertain so we use them combined: a C o A = 2. 



A. 2. 4. Double spiral-in 

If mass transfer is unstable when both stars are evolved 
(which can only happen if the mass ratio is close to unity) , 
we model the evolution as a common envelope in which 
the two cores spiral-in. The energy needed to expel the 
complete envelope is computed analogously to the case 



of a standard common envelope (Webbink 1984, see also 



Sect. A.2.2) 



Mi (Mi - M f ) TOi (m ; - m f ) 



A R 



Xr 



Mf TOf Mi TOi 

2 af 



2 Oi 



If the final separation is too small for the two cores to form 
a detached binary, the cores merge and we compute the 
fraction of the envelopes that is lost with the (practical) 
assumption that both stars lose the same fraction of mass, 
retaining /M, i.e 



M(l-/)Mi , TOi(l-/)TOi 



A R 



Xr 



/M frm Mi mi 
2 orlof 2 ft; 



where orlof is the separation at which one of the cores 
fills its Roche lobe. This is solved for /. 



A. 2. 3. Envelope ejection 



In the case of envelope ejection (Nelemans et al. 200C), we 
assume that the complete envelope is lost and that this 
mass loss reduces the angular momentum of the system 
linearly proportional to the mass loss, as first suggested 
for the general case of non-conservative mass transfer by 



Paczyhski & Ziolkowski (1967) 



Ji-J m =7^i^-, (A. 15) 



where Ji is the angular momentum of the pre-mass trans- 
fer binary and M tot is the total mass of the binary. 



